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In  designing  a  direct  methanol  fuel  cell  and  evaluating  the  appropriateness  of  new  materials  for  inclu¬ 
sion,  it  is  helpful  to  consider  the  impact  of  material  properties  on  the  performance  of  a  complete  system: 
to  some  degree,  poor  fuel  utilization  and  performance  losses  from  methanol  crossover  can  be  mitigated 
by  proper  system  design.  Simple  engineering  models  can  be  useful  tools  in  facilitating  this  type  of  system 
design.  In  this  paper,  an  analytical  model  is  developed  to  determine  the  oxygen  concentration  profile  in 
the  cathode  backing  layer  and  flow  channel  along  a  one-dimensional  cross-section  of  the  fuel  cell.  An 
existing  analytical  model  is  then  used  to  determine  the  methanol  concentration  profile  in  the  anode 
backing  layer  and  membrane  and  the  methanol  crossover  current  along  the  same  cross-section.  Applying 
a  fixed  cell  potential  and  using  Tafel  kinetics  to  describe  the  charge-transfer  reactions  at  the  anode  and 
the  cathode,  the  local  current  density  and  the  rates  of  methanol  and  oxygen  consumption  are  deter¬ 
mined.  The  process  repeats  for  a  fixed  number  of  cross-sections  down  the  entire  length  of  the  flow 
channels,  and  the  average  current  density  and  stoichiometric  ratios  are  calculated  at  the  end.  The  model 
is  applied  to  examine  the  effects  of  new  low-crossover  membranes  and  to  suggest  new  design  param¬ 
eters  for  those  membranes.  Also,  an  analysis  is  presented  in  which  the  tradeoff  between  stack  power  and 
the  size  of  system  components  is  examined  for  a  range  of  methanol  and  oxygen  flow  rates. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  an  electrochemical 
energy  conversion  device  that  operates  by  converting  methanol 
and  water  to  carbon  dioxide  and  protons  at  the  anode, 
CH3OH  +  H20  ^  CO2  +  6H+  +  6e  ,  and  by  reducing  oxygen  at  the 
cathode,  O2  +  4H+  +  4e”  2H20.  The  primary  advantage  of  the 

DMFC  is  that  it  can  be  operated  with  a  liquid  anode  feed  and 
therefore  does  not  require  a  reformer  or  a  bulky  gas  storage  tank. 
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Also,  because  the  oxidation  reaction  produces  six  electrons  per 
methanol  molecule,  the  DMFC  delivers  a  high  specific  energy.  The 
simplicity  and  compactness  of  the  fuel  delivery  system  makes  the 
DMFC  attractive  for  portable  applications. 

A  paramount  problem  in  the  implementation  of  the  DMFC  is 
methanol  crossing  over  to  the  cathode,  where  it  is  oxidized, 
resulting  in  depolarization  of  the  positive  electrode  and  wasted 
fuel.  The  origin  of  crossover  lies  in  the  fact  that  methanol  is 
completely  soluble  in  water  and  that  water  is  needed  to  swell  the 
proton-exchange  membrane  and  to  impart  protonic  conductivity 
across  the  membrane.  One  way  to  mitigate  this  problem  is  to  create 
a  low-crossover  membrane  to  replace  standard  perfluorosulfonic 
acid  membranes  such  as  Nafion®.  Increasing  the  efficiency  of  the 
reactions  at  the  anode  and  the  cathode  through  the  design  of  the 
catalyst  layers  can  also  minimize  the  harmful  effects  of  crossover. 
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Finally,  it  is  important  to  optimize  the  design  of  the  flow  channels 
and  operating  parameters  such  as  the  methanol  and  oxygen  feed 
concentrations  and  the  temperature. 

During  the  past  two  decades,  numerous  approaches  to  modeling 
the  DMFC  have  been  taken.  Most  of  these  models  focus  on  specific 
aspects  of  the  fuel  cell,  such  as  polarization,  multiphase  transport, 
or  methanol  crossover.  In  the  majority  of  the  models,  fuel  utiliza¬ 
tion,  system  efficiency  and  other  system-level  engineering  prob¬ 
lems,  if  considered  at  all,  is  only  given  brief  discussion.  Flowever, 
there  are  still  many  DMFC  models  in  the  literature  that  focus 
explicitly  on  issues  such  as  fuel  utilization  [1—7],  system  efficiency 
[8,9],  optimization  of  methanol  and  oxygen  flow  throughout  the 
fuel  cell  stack  [10-14],  exergy  and  thermal  analyses  [14-17],  and 
system  size  and  cost  [11,16,18],  What  most  of  these  models  lack  are 
simple  connections  among  materials  properties,  operating 
parameters,  and  system  constraints. 

This  work  seeks  to  make  those  connections  by  incorporating  as 
many  analytic  expressions  as  possible.  Specifically,  the  goal  here  is 
to  examine  the  effect  of  materials  properties,  especially  properties 
relating  to  the  proton-exchange  membrane,  on  the  net  power 
output  and  total  size  of  the  fuel  cell  system.  System  size  is  a  critical 
factor  for  the  cost  and  portability  for  small  DMFCs.  Running  at  high 
utilization  rates  (low  flow  rates)  will  lower  the  size  of  the  system 
components  such  as  the  methanol  tank  and  the  air  blower  but  will 
reduce  performance  and  require  a  larger  fuel  cell  stack  to  achieve 
a  certain  power  output.  Conversely,  running  at  low  utilization  rates 
(high  flow  rates)  will  increase  power  output  and  reduce  the  stack 
size  but  will  require  a  larger  methanol  tank  and  air  blower.  To 
perform  this  analysis,  we  designed  a  model  that  incorporates  as 
many  analytic  solutions  as  possible  for  simple  use  and  fast  run 
times  in  a  MATLAB  simulation  environment. 

To  create  the  desired  model,  we  seek  to  derive  analytical  solu¬ 
tions  for  the  methanol  and  oxygen  concentration  profiles  along 
a  single  one-dimensional  cross-section  of  the  fuel  cell.  The  solu¬ 
tions  for  the  methanol  concentration  in  the  anode  backing  layer 
and  membrane  taken  from  our  previous  publication  [11],  Similar 
equations  are  derived  in  this  work  for  the  oxygen  concentration  in 
the  cathode  backing  layer  and  flow  channel.  We  then  apply  Tafel 
kinetics  with  a  fixed  cell  potential  to  determine  the  local  current 
density  and  the  consumption  of  methanol  and  oxygen  reactants 
along  a  single  cross-section  of  the  cell.  Using  a  forward  difference 
approach,  the  methanol  and  oxygen  concentrations  at  the  next 
point  in  the  flow  channel  are  calculated.  The  model  is  stepped 
down  the  flow  channels  in  this  manner,  and  the  average  current 
density  is  calculated  by  dividing  the  sum  of  all  the  cross-sectional 
currents  by  the  number  of  steps. 

In  the  first  part  of  the  analysis,  we  simulate  polarization  data  for 
a  cell  operating  at  fixed  inlet  concentrations  and  flow  rates  of 
methanol  and  oxygen.  In  order  to  validate  the  model,  we  compare  the 
model’s  predictions  with  experimental  data  for  the  performance  of 
DMFCs  running  at  low  utilization  rates  of  both  oxygen  and  air  at  65  °C 
and  90  °C.  This  data  provides  a  base  case  for  setting  the  standard 
kinetic  parameters  in  the  model,  such  as  the  exchange  current 
densities  and  the  symmetry  coefficients.  We  then  apply  the  model  to 
anticipate  the  effects  of  using  new  low-crossover  membranes  from 
Zhu  et  al.  on  fuel  cell  performance,  and  we  offer  suggestions  on  how 
to  design  and  operate  fuel  cells  that  use  the  new  membranes. 

In  the  second  part  of  the  analysis,  we  consider  the  system-level 
tradeoff  between  the  power  output  of  the  cell  and  the  size  of  the 
methanol  storage  tank  and  air  blower.  We  set  a  desired  current 
output  for  the  cell  (50  mA  cm  2  or  300  mA  cm  2)  and  calculate  the 
power  output  at  different  oxygen  and  methanol  flow  rates  and 
methanol  concentrations.  We  also  consider  the  system-level 
benefits  of  using  the  new  membranes.  The  size  of  the  methanol 
tank  and  the  air  blower  should  scale  linearly  with  the  flow  rates 


while  the  marginal  increase  in  power  output  will  diminish  as  the 
flow  rates  increase.  Therefore,  given  a  set  of  performance 
constraints,  we  expect  to  find  flow  rates  that  minimize  the  total  size 
of  the  fuel  stack  in  combination  with  the  other  system  components. 

2.  Model  development 

2.1.  Assumptions 

In  Section  2.1,  we  develop  analytic  solutions  for  the  oxygen 
concentration  in  the  cathode  backing  layer  and  across  the  thickness 
of  the  cathode  flow  channel.  For  the  analytic  solutions  of  the 
methanol  concentration  profile  in  the  anode  backing  layer  and  in 
the  membrane,  we  again  refer  the  reader  to  our  previous  publica¬ 
tion  [4],  The  notation  from  that  publication  is  carried  into  this  work. 
Section  2.2  shows  how  the  local  current  density  is  found  using  the 
Tafel  equations  and  assuming  a  constant  cell  voltage.  In  Section  2.3, 
we  use  the  local  current  density  to  calculate  the  local  methanol  and 
oxygen  consumption  in  the  flow  channels  and  then  apply  a  forward 
difference  method  to  determine  the  methanol  and  oxygen 
concentration  profiles  down  the  entire  length  of  the  channels.  At 
the  end,  the  stoichiometric  coefficients  are  calculated  using  the 
average  current  density  in  the  fuel  cell. 

In  the  development  of  the  model  that  follows,  we  assume  that  all 
fluxes  and  current  densities  are  defined  per  unit  of  superficial  area, 
and  transport  properties  in  the  backing  layers  and  membrane  are 
effective  transport  properties,  modified  from  bulk  values  by 
appropriate  corrections  for  volume  fraction  and  tortuosity.  The  term 
“current”  will  always  refer  to  the  specific  current  density  unless 
otherwise  noted.  We  also  neglect  convection  across  the  thickness  of 
the  cell.  The  methanol  flux  across  the  cell  is  dominated  by  diffusion 
in  the  anode  backing  layer  and  is  a  combination  of  diffusion  and 
electroosmotic  drag  in  the  membrane,  while  the  oxygen  flux  is  due 
solely  to  convection  in  the  flow  channel  and  molecular  diffusion  in 
the  cathode  backing  layer.  The  methanol  concentration  is  driven  to 
zero  at  the  boundary  between  the  membrane  and  the  cathode 
backing  layer,  whereas  the  oxygen  partial  pressure  maintains 
a  finite  value  at  that  same  boundary  (Fig.  1 ). 


diffusion  diffusion 

only  and  drag  no  methanol 


N  / 

thin  catalyst  layers 


Fig.  1.  Cross-section  of  the  cell  including  a  qualitative  rendering  of  the  methanol  and 
oxygen  concentration  profiles. 
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The  most  critical  assumption  in  this  model  is  that  any  effects 
stemming  from  water  transport  or  carbon  dioxide  gas  removal  are 
ignored  or  incorporated  into  effective  diffusion  coefficients.  This 
assumption  allows  for  analytic  solutions  for  the  oxygen  concen¬ 
tration  to  be  developed,  but  it  also  restricts  the  validity  of  the  model 
to  low-power  applications.  In  systems  using  high  current  density 
and  high  methanol  concentration  (above  2  M),  or  in  any  scenario 
where  water  flux  in  the  cathode  and  carbon  dioxide  generation  in 
the  anode  and  cathode  are  substantial,  the  model’s  ability  to 
provide  accurate  performance  predictions  is  limited.  However,  in 
applications  such  as  portable  communications,  where  system  effi¬ 
ciency  and  system  size  are  paramount,  the  fuel  cell  is  often  oper¬ 
ated  at  low  current  densities.  These  systems  need  to  run  lean  on 
fuel  in  order  to  carry  the  methanol  oxidation  reaction  to  comple¬ 
tion,  thereby  avoiding  wasted  fuel  and  the  production  of 
intermediaries.  In  these  low-power  systems,  the  model  should 
have  a  wide  range  of  applicability. 

Another  important  simplification  is  the  assumption  of  Tafel 
kinetics  for  both  methanol  oxidation  and  oxygen  reduction.  Both 
reactions  are  multi-step,  and  Tafel  kinetics  only  applies  to  the  rate¬ 
determining  step  in  each  case.  However,  in  the  cases  described 
here,  the  exchange  current  densities  for  both  reactions,  ~10  5  for 
methanol  oxidation  and  ~  10  7  for  oxygen  reduction,  are  at  least 
two  orders  of  magnitude  below  the  cell  current  density,  which 
ranges  from  about  1  mA  cm-2  to  1  A  cm-2  in  the  experimental  data 
used  to  verify  the  model.  Therefore,  the  reaction  is  forward-biased 
to  such  an  extent  that  we  can  assume  the  intermediate  steps 
proceed  rapidly  and  apply  the  Tafel  equation  to  the  rate¬ 
determining  step.  As  long  as  the  current  density  is  not  so  high 
that  adsorption  of  the  reactants  on  the  catalyst  sites  becomes  rate¬ 
determining,  Tafel  kinetics  should  be  sufficient  to  describe  the 
entirety  of  the  electrode  reactions.  See,  for  example,  an  analysis  of 
methanol  reaction  kinetics  under  relevant  conditions  in  Ref.  [26]. 

The  methanol  solution  is  assumed  to  be  an  ideal  incompressible 
liquid  in  an  isothermal  environment  such  that  any  pressure  or  local 
temperature  effects  can  also  be  ignored.  The  air  in  the  cathode  is 
assumed  to  be  an  ideal  gas,  and  only  water  in  the  vapor  phase  is 
considered  explicitly  at  the  cathode.  Much  work  has  been  done,  most 
notably  by  Scott  et  al.  [19]  and  by  Wang  and  Wang  [20],  to  show  that 
two-phase  effects  at  the  anode  can  be  important  even  in  liquid-feed 
DMFCs,  but  this  model  does  not  consider  those  effects.  Instead,  we 
assume  that  the  behavior  of  the  backing  layers  can  be  characterized 
simply  by  an  effective  diffusion  coefficient  and  by  their  thickness. 


^H2o,sat  =  exp  11.6832 


3816.44  1 
(T  -  46.13)J  ’ 


0) 


then  the  oxygen  partial  pressure  will  be  20.946%  of  the  remaining 
total  pressure  for  an  air  feed.  A  study  by  Weber  and  Newman  [22] 
confirms  that  when  the  cathode  backing  layer  in  a  PEM  fuel  cell  has 
a  high  porosity  and  a  relatively  small  pressure  gradient  across  it, 
convection  and  Knudsen  diffusion  can  be  safely  ignored.  Therefore, 
in  order  to  obtain  an  analytic  solution  of  the  oxygen  concentration 
profile  across  the  backing  layer,  we  use  Stefan— Maxwell  diffusion 
to  fully  describe  the  transport  of  oxygen  in  that  region 


Again  note  that  D^_N2  anc*  ^o^-h2o  are  effective  diffusion  coef¬ 
ficients  for  the  components  in  the  porous  backing  layer.  The  flux  of 
water  and  nitrogen  across  the  backing  layer  is  negligible,  allowing 
us  to  ignore  the  last  two  terms.  Substituting  partial  pressures  for 
the  mole  fractions  and  N^ck  =  (i  +  ix)/4F  for  the  oxygen  flux,  we 
find  an  expression  for  the  change  in  absolute  oxygen  pressure 
across  the  backing  layer 


_  RT  ( _  Pn2  _  Ph2o  \  f i  +  ix\ 
0*  “  f total  ^  bgf  Nj  D*_Hj0J  v  4 F  )’ 


(3) 


where  Ptota]  =  P^ck  +  Pn2  +  Ph2o  with  PN2  and  PH2o  constant.  PH2q 
is  always  set  to  the  saturated  value  because  there  is  assumed  to  be 
some  liquid  water  in  the  cathode  backing  layer  at  all  times,  and  PN2 
is  either  79.054%  of  the  remaining  total  pressure  for  an  air  feed  or 
zero  for  an  oxygen  feed. 

Since  P^ack  is  included  in  Ptotai.  it  is  now  on  both  sides  of  this 
equation,  which  means  that  the  pressure  profile  in  the  backing 
layer  is  nonlinear.  The  appropriate  boundary  conditions  are 
P^ck  =  Pq2S  atx  =  0andP&fk  =  P02C  atx  =  Lbackc,  where  P02  S 
is  the  oxygen  pressure  at  the  boundary  between  the  backing  layer 
and  the  flow  channel  and  P02,C  is  the  oxygen  pressure  at  the 
boundary  between  the  membrane  and  the  backing  layer.  After 
integrating  with  these  boundary  conditions,  we  find  the  expression 


Po2,c  = 


—  (Pn2  +PH2o)  + 


(Po2,s+Pn2+Ph2o)  -2(t  +  ix)  Up 


Pn2 

Off  n. 


Dfn, 


(4) 


2.2.  Oxygen  pressure  in  the  cathode 
2.2.1.  Cathode  backing  layer 

In  the  cathode,  we  assume  that  the  incoming  air  stream  is 
a  pure  mixture  of  oxygen,  nitrogen  (if  air  is  used  instead  of  pure 
oxygen),  and  water  vapor  (if  humidified).  If  we  assume  that 
cathode  feed  stream  obeys  the  ideal  gas  law,  then  c^  =  Pq2  /RT, 
where  Pq2  =  yo2  Ptotai  >s  the  absolute  oxygen  pressure.  We  again 
emphasize  that  the  effects  of  any  liquid  water  present  in  the 
backing  layer  or  the  flow  channel  are  incorporated  into  the 
effective  diffusion  coefficients.  If  the  feed  stream  is  humidified 
and  if  the  saturation  vapor  pressure  of  water  is  solely  a  function 
of  temperature  [21] 


2.2.2.  Oxygen  transport  from  cathode  flow  channel  to  backing  layer 
A  crucial  difference  between  the  treatment  of  the  anode  and  the 
cathode  is  the  need  to  consider  the  movement  of  oxygen  gas  from 
the  bulk  of  the  flow  channel  to  the  surface  of  the  cathode  backing 
layer.  The  change  in  oxygen  partial  pressure  across  this  distance 
takes  the  same  form  as  Eq.  (3) 


3Pgfv  RT  /  PNa  Ph2o  \(i+ix\ 

0#  .Ptotai  V  Do2-N2  D02-H2O/  V  4F  )  ’ 


(5) 


where  Pg)nv  is  the  absolute  oxygen  pressure  in  the  convective 
region  between  the  bulk  flow  channel  and  the  backing 
layer  surface.  Note  that  the  diffusion  coefficients  lack  the  eff 
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subscript  because  we  are  outside  the  porous  region  of  the  backing 
layer. 

As  with  P*gck,  P“nv  is  included  in  Ptotai  and  is  on  both  sides  of  Eq. 
(5),  so  the  result  is  a  nonlinear  pressure  profile.  The  appropriate 
boundary  conditions  are  Pgnv  =  Pq2/  at  x  —  0  and  Pq”1™  — 
po  2  s  atx  =  dh/Shp,  i.e.  the  backing  layer  surface,  where  df,  is  the 
hydraulic  diameter  of  the  channel  and  Sh?  =  2.71  is  the  Sherwood 
number  for  a  square  channel  [23],  Pq2  j  is  the  oxygen  pressure  in  the 
bulk  flow  channel  and  is  determined  according  to  the  solution 
method  given  in  Section  2.4.2.  Integrating  Eq.  (5)  with  these 
boundary  conditions  results  in  an  expression  for  the  oxygen  pres¬ 
sure  at  the  surface  of  the  cathode  baking  layer  that  is  almost 
identical  in  form  to  the  expression  for  P02C. 


cathodic  activation  overpotential.  Because  of  the  complexity  of 
oxygen  diffusion  in  the  catalyst  layer,  which  has  relatively  low 
porosity  and  high  water  saturation,  we  choose  to  use  an  effective 
diffusion  coefficient  of  oxygen  d  based  on  the  idea  that  the 
catalyst  layer  is  fully  flooded  and  that  the  diffusion  of  oxygen  only 
occurs  in  the  regions  containing  Nafion  [24] 

<.d  =  D0  2, Nafion  [(1  ecl)  ^Nation]  •  (9) 

where 

Nafion  =  5.23  x  10-6exp[(23849/R)*((l/333.15)  _  (1/T))]cm2 
s-1  [25],  t'd  is  the  porosity  of  the  catalyst  layer,  and  <[)Naflon  is  the 
volume  fraction  of  Nafion  in  the  catalyst  layer.  Making  the  gov- 


P02,s 


-  (PN2  +  Ph2o)  3-  \j  (p0j.f  ■  PNj  •  P|l2o)  — 2(i  +  ix)  (jf  ^  (dqN2n 


Ph2Q  \ 

Dqj-HjO/  ShF 


(6) 


We  then  substitute  Pq2 s  directly  into  Eq.  ...  to  obtain  the  erning  equation  dimensionless  leads  us  to  define  the  Thiele 
following  expression  for  Pq2  C  terms  of  Pq2j  modulus  as  follows  [1  ]: 


C  =  —  (f*N2  +  ph2o)  +  ^  (Po2f  +  Pn2  +  Ph2o)  -2(i  +  ix) 


Pn2 

D02-n2 


Ph2o 

D02-h 


Pn2 


f  PH20 
Df_„ 


The  resulting  expression  cannot  be  simplified  and  must  be 
substituted  in  this  form  when  we  solve  for  the  desired  current  i. 
Aside  from  i,  the  other  unknown  in  Eq.  (7)  is  the  oxygen  pressure  in 
the  bulk  flow  channel  P02/.  The  solution  methods  for  i  and  P0lf  are 
described  in  Sections  2.2  and  2.3.2. 


2.2.3.  Oxygen  transport  in  the  cathode  catalyst  layer 

While  maintaining  the  assumption  of  thin  catalyst  layers,  we 
will  choose  to  define  an  effectiveness  factor  to  account  for  the  mass 
transport  of  oxygen  in  the  cathode  catalyst  layer  because  the 
oxygen  reduction  reaction  is  first-order.  The  methanol  oxidation 
reaction  is  still  assumed  to  be  zero-order,  so  there  is  no  need  to 
define  an  effectiveness  factor  there.  The  effectiveness  factor  is 
a  dimensionless  quantity  that  relates  the  actual  rate  of  reaction  to 
the  rate  if  the  entire  interior  of  the  catalyst  layer  was  exposed  to  the 
conditions  that  exist  at  the  catalyst  layer/backing  layer  boundary.  In 
effect,  it  is  a  measure  of  how  far  into  the  catalyst  layer  the  oxygen 
penetrates  and  is  a  simpler  and  more  convenient  method  of 
accounting  for  mass  transport  in  the  catalyst  layer  than  applying 
a  full  model  containing  terms  for  molecular  diffusion,  Knudsen 
diffusion  (particularly  important  in  catalyst  layers  and  micro- 
porous  layers),  and  convection. 

Assuming  Tafel  kinetics,  at  steady-state  the  governing  equation 
for  the  material  balance  of  oxygen  in  the  catalyst  layer  is  [1] 


d2c(x)  if?  RT  | 
dx2  4RclDeffd  P02,ceXP| 


( ^cathode  )j  c(x)  —  0- 


(8) 


where  if?  is  the  exchange  current  density  in  A  cm-2,  Lc i  is  the  catalyst 
layer  thickness,  Pq2  C  is  the  oxygen  pressure  at  the  catalyst  layer 
surface,  ac  is  the  cathodic  transfer  coefficient,  and  VCathode  >s  the 


no) 

Subsequently,  the  effectiveness  factor  is  a  function  of  the  Thiele 
modulus:  p  =  (tan/i  8)16.  For  an  effectiveness  factor  close  to  unity, 
internal  mass  transport  in  the  catalyst  layer  is  minimal,  and  the 
reactant  penetrates  fully  into  the  catalyst  layer.  If  the  effectiveness 
factor  is  closer  to  zero,  mass  transport  in  the  catalyst  layer  is 
a  significant  issue,  and  the  catalyst  layer  is  not  being  fully  utilized.  In 
DMFC  cathodes  the  effectiveness  factor  is  much  less  of  an  issue  than 
in  PEMFCs  due  to  lower  operating  current  densities.  Even  with  the 
relatively  thick  cathode  catalyst  layer  in  this  model  (Lc i  =  10  pm),  the 
effectiveness  factor  is  rarely  below  0.9.  However,  because  including 
it  does  not  add  noticeably  to  the  computation  time,  we  choose  to 
include  it  for  the  sake  of  completeness  and  accuracy. 

2.3.  Overpotentials  and  current  solution  method 

2.3.1.  Determination  of  overpotentials  and  cell  potential 

Now  that  the  concentration  profiles  for  methanol  and  oxygen 
have  been  developed  for  a  one-dimensional  cross-section  of  the  cell 
in  terms  of  the  desired  current  i,  all  that  remains  is  to  solve  for  that 
current.  As  stated  in  the  introduction,  our  strategy  is  to  assume  the 
cell  potential  remains  constant  and  numerically  solve  for  the  local 
current  along  each  cross-section.  Given  the  high  in-plane  conduc¬ 
tivity  of  the  bipolar  plates,  any  voltage  difference  along  the  plates 
will  be  quickly  dissipated,  so  the  assumption  of  a  constant  cell 
potential  appears  reasonable.  The  cell  potential  is  defined  as  follows: 

'/cell  =  U6  -  L  -  '/anode  +  Vcathode,  (11) 

Kmem  ^  \Kcla  Kclc / 
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where  Kmem  is  the  membrane  conductivity  in  S  cm-1,  /ccia  and  /ccjc  are 
the  anode  and  cathode  catalyst  layer  conductivity  in  S  cm-1  and 
if  =  1.21  V  for  the  DMFC.  Implicit  in  this  equation  is  the  assumption 
that  the  proton  conductivity  is  uniform  throughout  the  thickness  of 
the  catalyst  layers  and  the  membrane.  While  each  is  tracked  sepa¬ 
rately,  and  preserves  the  ability  to  describe  restricted  proton  access 
in  the  catalyst  layers,  where  the  volume  fraction  of  ionomer  is  lower 
than  in  the  membrane  separator,  we  do  neglect  the  possibility  of 
variation  through  the  thickness  of  any  of  the  components. 
Furthermore,  the  factor  of  Vi  appears  in  front  of  the  catalyst  layer 
ohmic  terms  because  we  assume  that  the  protons  travel  on  average 
half  the  distance  across  each  catalyst  layer.  While  this  assumption  is 
not  true  in  a  situation  where  there  are  substantial  mass-transport 
limitations  in  the  catalyst  layers,  zero-order  kinetics  at  the  anode 
and  an  effectiveness  factor  that  is  almost  always  above  0.9  at  the 
cathode  support  the  notion  that  the  reactions  are  occurring  some¬ 
what  uniformly  across  the  catalyst  layers. 

As  mentioned  previously,  both  the  methanol  oxidation  and 
oxygen  reduction  reactions  are  assumed  to  follow  Tafel  kinetics, 
with  the  former  being  zero-order  and  the  later  being  first-order. 
The  assumption  of  zero-order  kinetics  has  been  shown  to  be 
accurate  to  very  low  concentrations  or  operation  very  near  the 
limiting  current,  and  is  tantamount  to  assuming  that  the  CO-sites 
on  the  anode  electrocatalyst  can  quickly  become  saturated  [26], 


Substituting  Eq.  (15)  into  Eq.  (16)  and  solving  for  Canode  produces 
a  relationship  between  canode  and  Cfeed 


Further  substitution  of  ilim,a  =  6F(D^fe0H_H20/Lbacka)cfeed 
results  the  following  expression  for  canode: 


_  ’/'Cfeed  ‘-feed 
Lanode  -  f(\+p)  ' 


(18) 


Eq.  (18)  shows  that  canode  is  dependent  on  Cfeed  at  the  corresponding 
location  in  the  flow  channel,  the  initial  feed  concentration  c®eed,  the 
dimensionless  ratio  v,  and  if  =  i}jma/i,  which  is  the  limiting  current 
at  the  anode  entrance  normalized  by  the  desired  current. 

Removing  f  from  the  denominator  and  substituting  the  above 
expression  for  canode  into  Eq.  (15)  results  in 


f  him, a  i  A^Cfeed  +  C°A 

<>+’)  j  (  1 

Noting  that  km,a/^ima  =  cfeed/cfeed’  we  arrive  at  the  desired 
relation  for  ix  in  terms  of  i. 


^anode  — 


.  /  V'f’Cfeed  +  Cjeed 

*  A  ^d(1+») 


(20) 


Note  the  presence  of  the  effectiveness  factor  in  the  cathodic 
overpotential.  Maintaining  zero-order  methanol  kinetics  means 
that  the  anode  overpotential  does  not  depend  on  the  methanol 
concentration  at  the  anode  catalyst  layer.  This  assumption  is  valid 
as  long  as  we  keep  the  methanol  stoichiometry  high  enough  to 
avoid  depleting  methanol  at  the  anode  catalyst  layer  and  incurring 
a  transition  to  first-order  methanol  kinetics  [26].  Substituting  Eqs. 
(12)  and  (13)  into  Eq.  (11 )  results  in  the  following  expression  for  the 
cell  potential: 


It  is  clear  that  we  must  use  a  numerical  method  to  solve  for  the 
current  in  Eq.  (14).  However,  in  order  to  find  the  desired  current,  we 
must  first  find  an  expression  for  the  crossover  current  ix  in  terms  of 
the  desired  current  i  and  other  known  variables. 


2.3.2.  Determination  of  crossover  current  and  total  current 

To  write  ix  in  terms  of  i  we  refer  to  our  previous  work  [11] 
concerning  the  methanol  concentration  in  the  anode  backing 
layer  and  in  the  membrane.  The  crossover  current  is  described  in 
terms  of  the  feed  concentration  Cfeed  and  the  concentration  at  the 
anode  catalyst  layer  canode  as  follows: 


The  crossover  current  can  also  be  described  using  the  dimen¬ 
sionless  parameter  v. 


ix  =  (W-0  (!+,)■  H6) 


Now  we  are  prepared  to  make  the  appropriate  substitutions  and 
solve  Eq.  (14)  for  i.  The  solution  is  performed  in  MATLAB  using  the 
root-finding  function  fzero  along  each  cross-section  of  the  cell.  To 
obtain  an  average  cell  current  the  sum  of  the  cross-sectional 
currents  is  divided  by  the  number  of  steps. 


2.4.  Consumption  of  methanol  and  oxygen  in  the  flow  channels 
2.4.1.  Anode  flow  channel 

Having  the  ability  to  find  the  local  current  density  as  a  function 
of  the  local  methanol  concentration  and  oxygen  pressure  in  the 
flow  channels,  we  now  must  determine  the  consumption  of 
methanol  and  oxygen  in  the  flow  channels  as  we  move  from  one 
step  to  the  next  (Fig.  2). 

Assuming  the  velocity  vm^oh  of  methanol  in  the  flow  channel  is 
constant,  the  total  flow  rate  of  methanol  in  the  anode  flow  channel 
is  given  by 


anode  cathode 

backing  backing 


Fig.  2.  Methanol  consumption  in  the  fuel  cell. 
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WMeOH  -  cfeed  Ahan  vMeOH •  (21) 

where  /\chan  =  wChanhchan.  The  change  in  feed  concentration 
between  each  step  as  we  move  down  the  flow  channel  is 

'•Cfeed  _  ]\fback  (  wsepLchan  A  _  /i'+ljA  (  wsepfchan  ^ 

-*r—  MeOH  Ute-teJ  “ — V6 r)  UchanvMe0HJ  •  (22) 

where  z'  =  z/Lc han  is  the  normalized  distance  down  the  channel,  and 
wsep  is  the  channel  separation  width,  that  is,  the  width  of  electrode 
perpendicular  to  the  channel  flow  that  is  fed  by  a  single  channel. 
Wseptchan  is  the  area  of  the  anode  that  receives  the  methanol 
flowing  through  the  channel.  The  effects  of  other  parallel  channels 
are  ignored. 

Noting  the  presence  of  the  crossover  current  term  ix,  we 
substitute  Eq.  (20)  to  obtain 


9Cfeed  V^Cfeed  -_c|ed.\  MepW\ 

9Z'  "  4^(1+,)  )\  iW  )■  ^ 

We  can  now  define  the  stoichiometric  ratio  for  the  methanol 


where  K  —  D/j/2/Re  represents  the  cross-sectional  area  to  flow.  Reis 
the  Reynolds  number, /is  a  friction  factor,  and  their  product  is  equal  to 
14.2296  for  a  square  channel  [23],  D/,  is  the  hydraulic  diameter  and  is 
equal  to  four  times  the  area  of  the  channel  divided  by  the  perimeter. 
/iair  is  the  dynamic  viscosity  of  the  air  mixture  and  is  a  function  of  the 
mole  fraction  of  all  three  components  in  the  mixture 


Pair 


im 


(28) 


where  N  is  the  total  number  of  species  in  the  mixture,  y,-  and  yj  are 
the  mole  fractions  of  species  i  and  j,  and  'I>,  is  a  sum  of  dimen¬ 
sionless  ratios.  For  an  in  depth  analysis  of  the  derivation  of  /zair,  we 
refer  the  reader  to  pp.  153-155  in  Ref.  [23],  At  65  °C,  fi3 ir  is  about 
1.6  x  10”5  Pa  s_1,  depending  on  the  values  of  the  mole  fractions. 

Substituting  for  k  and  noting  that  Nairj/P3irj  —  Nq2j/Pq2j.  we 
obtain  a  final  expression  for  the  pressure  drop  as  a  function  of  the 
oxygen  flow  rate,  oxygen  pressure,  and  known  parameters 


9Pair/  _  2/Re/tair  N0lJ  RT 
Df,  A:han  Po2f' 

Using  a  forward  difference  to  discretize  this  equation,  as  we  did 
with  Eq.  (25),  results  in 


sMeOH  M 


cfeed  Alum17 

i/6F  tchanwsep 


(24) 


and  derive  an  expression  for  the  change  in  the  feed  concentration 


aCfeed  =  f  1-  A  /Vfeed  ~  CfeecA 
9Z'  V SMeOH /  \  (l+y)  J 


(25) 


(30) 

Therefore  the  pressure  drop  at  step  k  +  1  is  determined  by  the 
ratio  of  the  oxygen  flow  rate  to  the  oxygen  pressure  at  step  k.  The 
change  in  oxygen  flow  rate  is  simply  a  function  of  the  total  current, 
and  that  expression  can  be  discretized  to  solve  for  N0lj(k). 


The  stoichiometric  ratio  is  the  rate  at  which  methanol  is  deliv¬ 
ered  to  the  inlet  of  the  cell  relative  to  the  amount  of  methanol  that  is 
consumed  in  the  preferred  electrochemical  reaction.  Hence,  it 
represents  the  inverse  of  fuel  utilization  or  fuel  efficiency.  In  order  to 
discretize  the  flow  channel  and  apply  the  one-dimensional  model 
from  Sections  2.1  to  2.3,  we  apply  a  first-order  forward  difference, 
J{x+h)  =J[x)  +  hfx,  where  h  —  Lchan/n  is  the  step  size  and  n  is  the  total 
number  of  steps.  Then  we  create  the  following  expression  for  Cfeed  at 
step  k  +  1  based  on  Cfeed  and  the  current  i  at  step  k. 

*«„<*+ „  -  (T)  (JJ  m 


2.4.2.  Cathode  flow  channel 

In  the  cathode  flow  channel,  where  the  oxygen  reactant  is  a  gas,  we 
must  consider  the  pressure  drop  of  the  air  stream  in  addition  to  the 
consumption  of  oxygen.  We  cannot  assume  that  the  velocity  of  the  air 
stream  is  constant  as  we  do  for  the  velocity  of  the  liquid  methanol/ 
water  solution  in  the  anode  flow  channel.  We  do  assume  that  the 
absolute  pressure  of  water  vapor  remains  constant  throughout  the 
flow  channel.  If  the  inlet  flow  is  dry,  then  the  water  vapor  pressure  is 
set  to  zero,  and  if  it  is  humidified,  then  the  water  vapor  pressure  is 
assumed  to  be  its  saturated  value  given  in  Eq.  ( 1 ).  This  simplification  is 
necessary  given  that  changes  in  water  content  are  not  treated  in  the 
backing  layer,  but  it  can  lead  to  significant  error  at  high  current  when 
the  water  generation  is  high.  Because  of  these  assumptions,  the 
pressure  drop  only  affects  the  nitrogen  and  oxygen  pressures.  In 
general  terms,  according  to  Darcy’s  Law,  the  pressure  drop  is 


dPair/  Jtajr  =  _Mair  Nair/ 
dZ  ~  K  K  Achan  Pair/- 


(27) 


%/('<+!)  =  ^  (31) 

Determining  the  local  oxygen  pressure  P^jik)  is  a  bit  more 
challenging.  The  expression  for  the  consumption  of  oxygen  in  the 
cathode  flow  channel  is  derived  in  a  similar  manner  as  for  methanol 
in  the  anode.  The  change  in  oxygen  pressure  due  to  consumption  is 


where  the  substitution  v0lj  —  v3irj  —  ( N0lf /Adl3n)(RT/P02j )  is 
made.  Substituting  Eq.  (20)  for  ix,  we  find  an  expression  for  the 
change  in  oxygen  pressure  that  is  similar  to  Eq.  (23) 


9Pq2/ 

az 


'Mfeed  ~  Cfeed^ 

V  J 


(33) 


Discretizing  Eq.  (33)  and  adding  it  to  Eq.  (30)  results  in  an 
expression  for  the  change  in  the  total  pressure  of  the  air  inlet 
stream  at  each  step 

p,V(M)  -  p -  (l^®)  Q 

Acha^  f2fRenair(k)\ 

WjW)  V  siz  )\  Dl  J 

(No2j(k)  RT  \ 

X{^r  p^m)' 


(34) 
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Because  the  total  pressure  is  the  sum  of  the  partial  pressures, 
with  the  water  vapor  pressure  being  constant,  we  can  relate  the 
nitrogen  and  oxygen  pressures  in  order  to  make  substitutions 
and  solve  for  each  partial  pressure  individually:  *V  =  Pn2/ 
(Wo2//%/)andPN2/=Po2/(NN2//N02/).  The  analysis  above 
applies  to  a  system  where  the  flow  channels  are  in  a  coflow 
configuration,  in  which  the  methanol  and  oxygen  flow  in  the  same 
direction  through  the  fuel  cell.  If  we  desire,  we  can  add  the  change 
in  oxygen  pressure  between  steps  instead  of  subtracting  it  and 
thereby  simulate  a  counterflow  configuration,  in  which  the  meth¬ 
anol  and  oxygen  flow  in  opposite  directions.  The  only  difficulty  in 
that  situation  is  the  need  to  define  the  oxygen  pressure  at  the 
channel  outlet  and  the  oxygen  flow  rate  at  the  inlet,  both  of  which 
change  as  we  move  down  the  flow  channel.  For  the  sake  of 
simplicity,  we  choose  to  use  a  coflow  configuration  in  the  analysis 
below.  The  oxygen  stoichiometric  ratio  remains  very  high  (always 
above  10)  in  the  experimental  data  used  for  validating  the  model, 
which  justifies  this  approach.  Furthermore,  the  flow  channels  have 
a  single-serpentine  design,  and  there  is  probably  a  large  amount  of 
under-rib  transport.  At  the  high  flow  rates  we  are  considering,  we 
expect  there  to  be  little  difference  between  coflow  and  counterflow 
in  the  current/potential  regimes  that  are  not  mass-transfer  limited. 

3.  Experimental  validation 

3.1.  Nafion  membranes  with  optimized  catalyst  layers 

The  purpose  of  this  model  is  to  quickly  and  easily  determine 
the  effect  of  materials  properties  and  operational  parameters, 
especially  the  membrane  permeability,  conductivity,  and  thickness 
and  the  methanol  and  oxygen  flow  rates,  on  the  performance  and 
the  required  size  of  the  DMFC  system.  In  order  to  fit  certain 
parameters  and  accurately  predict  the  performance  of  real  systems, 
we  need  to  validate  the  model  with  a  control  set  of  experimental 
data.  One  drawback  to  this  model  is  the  assumption  of  a  negligibly 
thin  catalyst  layer  and  simple  zero-order  (methanol)  or  first-order 
(oxygen)  Tafel  kinetics  for  the  complex  reactions  taking  place  at  the 
anode  and  cathode  catalyst  layers.  Hence,  the  anodic  and  cathodic 
exchange  current  densities  (i°  and  i°)  and  the  symmetry  coeffi¬ 
cients  ( oia  and  ac)  must  be  adjusted  to  fit  experimental  data.  After 
comparison  with  the  data,  we  found  that  the  addition  of  a  contact 
resistance  term  Rc  was  unnecessary.  The  base  values  that  we  use  for 
all  of  the  physical  properties  are  listed  in  Table  1  and  do  not  change 
in  the  following  sections  unless  otherwise  noted.  The  values  of  all 
the  kinetic  and  resistive  properties  are  listed  in  Table  2.  It  is 
important  to  note  that  we  assume  liquid  water  to  be  almost  as 
effective  as  a  solid  medium  at  preventing  gas  diffusion.  Therefore, 
we  account  for  the  effect  of  liquid  water  in  the  cathode  by  lowering 
the  porosity  of  the  cathode  considerably  relative  to  the  anode. 

As  an  initial  check,  we  used  a  base  case  scenario  (Nafion  115, 
250  seem  oxygen  flow  rate,  and  2  ml  min  1  methanol  flow  rate)  and 
varied  the  temperature  (65  °C  or  90  °C)  and  the  cathode  flow 
composition  (air  or  pure  oxygen).  The  experiments  were  performed 
with  a  1  M  MeOH  feed  at  the  anode  and  20  psi  of  backpressure  at 
the  cathode.  The  resulting  four  cases  for  a  single  cell  with  a  5  cm2 
MEA  are  shown  in  Fig.  1  below.  The  MEAs  were  fabricated  by 
spraying  the  respective  catalyst  layer  inks  onto  each  side  of  the 
Nafion  membrane  and  then  hot-pressing  the  carbon  cloth  backing 
layers  onto  each  side.  The  catalyst  loading  was  2.5  mg  cm-2  Pt/Ru 
on  the  anode  and  2.5  mg  cm-2  Pt  on  the  cathode  on  unsupported 
carbon  black.  This  validation  was  performed  to  determine  the 
accuracy  of  the  temperature-dependent  diffusion  coefficients, 
conductivities,  and  exchange  current  densities  as  well  as  to  verify 
our  assumption  of  a  constant  water  vapor  pressure  in  the  cathode. 
These  experiments  also  tested  the  validity  of  our  choice  to 


Table  1 

Physical  properties  for  experimental  validation. 


Parameter  (unit) 

Value 

Reference 

Do-,-h,o  (cm2  s"1) 

6.436  X  1CT6  *  r1 

823/p( in  atm)  O’Hayre  et  al.  [23] 

Aj2-n2  (cm's-1) 

4.2  x  10-7  *  T2-33  ■ 

l/p( in  atm)  O’Hayre  et  al.  [23] 

DMeOH-H20  (cm2  S  ’) 

Yaws  [30] 

Dmem  (cm2  S  *) 

4.9  x  10-6  *  e24161 

[1/333  -  i/T)  Kauranen  and  Skou  [31] 

ea  (dimensionless) 

0.7 

Wang  and  Wang  [20] 

ec  (dimensionless) 

0.35 

Ed  (dimensionless) 

0.1 

Assumed 

^Nafion  (dimensionless)  0.058 

Measured 

J-back .a  =  4ack ,c  (cm) 

0.02 

Assumed 

^mem  (cm) 

0.0128 

Assumed 

Ph2o  (g  cm-3) 

1.1603  -  5.371  x 

IQ-4  .  T  incropera  et  al.  [21  ] 

£h2o  (dimensionless) 

2.5 

Zawodzinski  et  al.  [32] 

^MeOH  (dimensionless)  £h2o*18/Ph2o 

Calculated 

wsep  (cm) 

0.2 

Ichan  (Cm) 

25 

Measured 

Achan  (cm  ) 

0.01 

Measured 

incorporate  the  effects  of  liquid  water  and  carbon  dioxide  into  an 
effective  diffusion  coefficient  rather  than  consider  those  compo¬ 
nents  explicitly. 

The  primary  discrepancy  between  the  model  and  the  experi¬ 
ments  exists  in  the  mass-transport  region,  especially  in  the  65  °C 
case  with  air.  We  believe  this  discrepancy  exists  mostly  because  our 
model  does  not  allow  for  varying  water  content  in  the  catalyst  layer, 
gas  diffusion  layer,  and  flow  channel  of  the  cathode.  High  current 
density  generally  corresponds  with  high  water  crossover  and 
generation  rates  and  lower  effective  oxygen  diffusion  coefficients, 
which  decreases  performance.  However,  accounting  for  these 
variations  would  have  eliminated  our  ability  to  derive  analytic 
expressions  for  the  gas  concentrations  across  the  thickness  of  the 
cell.  We  decided  that  a  modeling  scheme  which  can  match  the 
experimental  performance  over  the  current  and  potential  ranges  of 
interest  (generally  less  than  300  mA  cm-2)  with  a  simple  set  of 
assumptions  and  a  fast  run  time  (usually  less  than  1  min)  is  of 
greater  interest  to  the  engineering  problems  of  system  sizing  and 
stack  optimization. 

Experimental  results,  including  the  results  shown  in  Fig.  3,  often 
use  deliberately  high  flow  rates  in  order  to  maximize  performance. 
However,  in  a  real  system,  it  might  behoove  the  designer  to  take 
a  small  cut  in  performance  by  lowering  the  flow  rate  to  say, 
50  seem,  such  that  the  air  blower  and  the  compressor  consume  less 
power.  The  ability  of  this  model  to  accurately  predict  fuel  cell 
performance  for  a  variety  of  methanol  and  oxygen  flow  rates  will  be 
important  when  we  consider  optimizing  the  size  of  system 
components  in  Section  4. 

3.2.  Low-crossover  crosslinked  blend  membranes 

The  cases  above  all  employ  Nafion  membranes,  but  much  work 
has  been  done  in  recent  years  to  develop  new  membranes  with 
lower  methanol  crossover  than  Nafion  and  comparable  proton 
conductivity  and  mechanical  properties.  In  particular,  Zhu  et  al. 
have  reported  novel  crosslinked  membranes  based  on  a  blend  of 


Table  2 

Kinetic  and  resistive  properties  for  Nafion  membranes. 


Parameter  (unit)  Value 
i°  (A  cm-2)  3.5  x  1(T5 

i°  (A  cm-2)  3.0  x  10-7 

a,,  (dimensionless)  0.8 

ac  (dimensionless)  0.875 

Kmem  (S  Cm-1 )  0.1 1 1  *  e1268(1'303  -  'IT> 

Kelt,  (S  cm-1)  0.001  *  e1268'1'298  - 

kc1c  (S  cm-1)  0.0026  *  e1268™298  “  ™ 


Reference 
Fitted  to  data 
Fitted  to  data 
Fitted  to  data 
Fitted  to  data 
O’Hayre  et  al.  [23] 

Havranek  and  Wippermann  [33] 
Saab  et  al.  [34] 


■etal./  j 


Fig.  3.  Comparison  of  simulations  to  experimental  data  at  65  °C  and  90  °C.  The  base 


ie)a 


that  showed  a  dramatic  improvement  in  performance  compared  to 
Nation  [27],  They  kindly  made  that  data  available  in  tabular  form 
along  with  physical  property  data  for  our  use  in  Figs.  4  and  5.  This 
section  describes  a  case  study  in  which  we  compared  our  model 
with  the  data  for  their  C-10  membrane.  These  tests  were  run  at 
65  °C  and  1  M  MeOH  using  a  humidified  oxygen  feed  at  a  200  seem 


Fig.  4.  Comparison  of  experiments  and  simulations  for  the  C-10  blend  membrane  and 
Nation  115  membrane  at  65  °C.  This  scenario  uses  a  2.5  ml  min-1  methanol  flow  rate 


etal.  [13], 


methanol  crossover.  There  are  several  experiments  in  the  literature, 
most  notably  from  Ge  and  Liu  [28],  running  under  similar  condi¬ 
tions  that  show  much  better  performance  from  Nafion  membranes 
than  what  is  shown  in  Fig.  4.  We  believe  the  discrepancy  at  higher 
currents  is  due  to  mass-transport  problems  in  the  experiment, 
specifically  water  management,  for  which  our  model  does  not 
account.  Also,  a  different  construction  of  the  catalyst  layers,  as  was 
done  with  the  experiments  shown  in  Fig.  3,  might  lead  to  better 
performance  for  both  the  Nafion  and  the  blend  membranes.  Finally, 
the  discrepancies  at  low  current  could  be  the  result  of  the  methanol 
diffusion  coefficient  in  the  blend  membranes  being  higher  during 
the  experiment  than  what  was  measured  prior  to  the  experiment. 

One  advantage  of  the  blend  membranes  is  that  they  can  be  made 
thinner  without  significantly  increasing  the  methanol  crossover, 
which  can  more  than  compensate  for  the  decrease  in  specific 
conductivity  when  compared  to  Nafion.  Fig.  5  simulates  the 
performance  of  the  C-10  membrane  reported  by  Zhu  et  al.  [27]  for 
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a  variety  of  thicknesses  under  the  conditions  given  in  Section  3. 
Fig.  6  does  the  same  with  Nation  membranes.  While  the  thinner 
Nation  membranes  suffer  from  noticeable  performance  declines  at 
low  current  due  to  crossover,  the  blend  membranes  show  no  drop 
in  performance.  As  we  mentioned  above,  the  measured  methanol 
diffusion  coefficient  might  be  lower  than  the  actual  value,  but  the 
model  still  clearly  shows  that  the  blend  membranes  should  be 
made  as  thin  as  manufacturing  tolerances  allow.  The  reduction  in 
ohmic  resistance  more  than  compensates  for  the  small  increase  in 
methanol  crossover. 

Aside  from  improved  performance,  the  lower  methanol  cross¬ 
over  rates  of  the  new  membranes  results  in  improved  fuel  effi¬ 
ciency.  If  we  want  to  calculate  the  maximum  fuel  efficiency  of  the 
fuel  cell  we  can  perform  an  iterative  reduction  of  the  methanol 
stoichiometric  ratio  until  the  methanol  concentration  at  the  anode 
reaches  the  first-order  transition  value.  Also  of  importance  when  it 
comes  to  optimizing  the  DMFC  system  with  new  membranes  are 
the  properties  of  the  backing  layers.  The  thickness  and  porosity  of 
the  backing  layers  must  strike  a  balance  between  allowing  reac¬ 
tants  to  diffuse  to  the  catalyst  layers  and  providing  an  effective 
barrier,  both  physically  and  chemically,  between  the  reaction  zones 
and  the  flow  channels.  The  anode  backing  layer  is  of  particular 
importance  because  it  must  allow  enough  methanol  to  reach  the 
anode  catalyst  layer  so  as  to  promote  the  desired  reaction  while  not 
allowing  so  much  methanol  so  as  to  promote  crossover.  Our 
previous  work  [4]  addresses  in  detail  the  issues  of  fuel  efficiency 
and  backing  layer  properties,  so  in  order  to  avoid  repetition,  this 
paper  will  not  discuss  these  issues  at  length. 

4.  The  effect  of  methanol  and  oxygen  flow  rates  on  system 
size  and  performance 

Because  portability  is  of  primary  importance  in  many  applica¬ 
tions  where  DMFCs  might  be  used,  we  want  to  use  this  model  to 
compare  the  power  consumption  of  the  air  blower  and  compressor 
to  the  stack  power  output  for  a  range  of  methanol  and  air  flow  rates. 
The  less  methanol  and  oxygen  that  are  used  and  the  more  slowly 
those  reactants  are  passed  through  the  flow  channels,  the  smaller 
we  can  make  the  supporting  system  components.  Ideally,  there 
would  be  no  methanol  crossover,  and  all  of  the  reactants  in  the 
system  would  either  be  consumed  in  the  desired  reaction  or  exit 
the  flow  channels,  perhaps  being  recirculated  through  the  system 
to  be  consumed  later.  In  this  case,  the  size  of  the  supporting  system 
components  would  be  solely  dependent  on  the  demand  for  current. 


The  parameters  and  operating  conditions  are  identical  to  what  is  used  for  the  Nafion 
115  membrane  in  Fig.  3. 


However,  the  existence  of  crossover  and  the  potential  for  wasted 
fuel  means  that  the  size  of  the  tank  and  the  blower  also  depends  on 
the  fuel  efficiency. 

If  reducing  the  size  of  the  system  components  were  our  only 
concern,  then  we  would  like  to  run  the  cell  as  lean  as  possible  with 
a  very  low  stoichiometric  ratio.  However,  running  a  lean  feed  will 
result  in  performance  losses  compared  to  running  at  a  high  stoi¬ 
chiometric  ratio.  We  are  especially  concerned  with  avoiding  the 
transition  from  zero-order  to  first-order  kinetics  for  the  methanol 
oxidation  reaction.  Rather  than  model  this  complicated  transition, 
which  would  require  a  substantial  increase  in  overpotential  at  the 
anode  and  a  corresponding  loss  in  performance  [26],  we  ensure 
that  the  optimized  system  remains  above  any  such  transition  by 
requiring  the  feed  concentration  to  remain  above  a  minimum  value 
throughout  the  flow  channel.  In  doing  so,  we  impose  a  constraint 
on  the  minimum  allowable  stoichiometric  ratio  at  the  anode.  At  the 
cathode,  our  only  constraint  on  the  stoichiometric  ratio  is  that  it 
remains  above  unity. 

It  is  important  to  note  that  the  analysis  in  this  section  does  not 
constitute  an  optimization  study  in  the  true  spirit  of  the  term.  We 
iterate  through  a  defined  set  of  flow  rates  for  a  few  specific 
scenarios  instead  of  using  an  optimization  function  to  explore  the 
entire  parameter  space.  A  more  thorough  optimization  analysis 
could  be  done  with  this  model,  but  we  choose  not  to  do  so  here  for 
the  sake  of  brevity.  Hence,  the  more  accurate  term  to  describe  what 
is  done  here  is  a  sensitivity  analysis.  The  goal  of  this  section  is  to 
illustrate  how  the  model  can  be  applied  across  a  wide  range  of 
a  given  parameter  space  and  how  it  can  provide  guidelines  for 
certain  parameters. 

For  a  single  fuel  cell,  the  size  of  the  fuel  tank  Vtan  k  scales  linearly 
with  the  flow  rate  at  the  anode  inlet  ZVj^Q^,  the  time  of  operation  t 
and  the  molar  volume  of  the  methanol  and  water  feed  stream 
^feed  :  ^tan  k  =  ^MeOH^feedT  where  NM£0H  =  c[?eed  FMeOH  as 
in  Eq.  (21 ).  We  have  already  seen  in  Fig.  6  how  the  performance  of 
the  fuel  cell  is  largely  independent  of  the  anode  flow  rate  due  to 
zero-order  kinetics,  so  in  that  sense,  we  want  to  run  at  as  low  a  flow 
rate  as  possible  to  avoid  crossover  and  wasted  fuel.  However,  there 
must  be  a  point  at  which  the  flow  rate  is  not  high  enough  to  sustain 
zero-order  kinetics.  Using  the  model  to  predict  this  minimum  flow 
rate  would  be  useful  in  determining  parameters  for  fuel  cell  design 
and  operation. 

The  scenario  presented  in  Fig.  7  above  compares  the  power 
output  of  a  fuel  cell  operating  at  300  mA  cm-2  current  using  Nafion 
115  to  fuel  cells  using  blend  membranes  of  the  same  thicknesses 
(represented  by  the  blue  lines)  as  in  Fig.  5.  The  operating  conditions 


operating  conditions  are  identical  to  what  is  used  in  Fig.  3. 
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and  the  cell  design  are  the  same  for  the  Nafion  and  blend 
membranes.  The  blend  membranes  are  very  effective  at  preventing 
methanol  crossover,  so  the  model  tells  us  that  we  should  manu¬ 
facture  them  as  thin  as  possible  in  order  to  lower  the  ohmic 
resistance,  especially  at  high  current.  Of  course,  the  limits  of  the 
manufacturing  process  and  the  desire  to  maintain  good  mechanical 
properties  and  durability  will  limit  our  ability  to  reduce  the 
membrane  thickness. 

In  this  case,  the  minimum  methanol  flow  rate  is 
4  x  10  6  mol  s-1,  which  corresponds  to  0.25  ml  min-1  and  an 
average  methanol  stoichiometric  ratio  of  1.6.  Since  the  model 
cannot  handle  the  transition  to  first-order  kinetics,  the  three  points 
below  that  transition  are  generated  by  lowering  the  desired  current 
until  the  methanol  concentration  remains  above  the  transition 
value  all  the  way  to  the  anode  exit.  Hence,  the  power  output  of  the 
cell  must  decrease  either  by  incurring  a  drop  in  cell  voltage  due  to 
the  transition  to  first-order  kinetics  or  by  lowering  the  current  until 
that  transition  is  avoided.  In  designing  a  system  under  our  given 
conditions,  we  would  want  to  keep  the  methanol  flow  rate  close  to 
this  value  but  safely  above  it  to  maximize  performance  and 
efficiency. 

The  situation  is  similar  when  we  consider  the  oxygen  flow  rate 
under  the  same  set  of  conditions.  Given  the  low  flow  rates  that  we 
are  using,  the  power  consumption  of  the  air  blower  and  air 
compressor  will  scale  linearly  with  the  flow  rate  in  the  cathode 
flow  channel.  Therefore,  as  with  the  anode  flow  rate,  we  would 
expect  at  to  see  a  point  at  which  the  extra  power  generated  by  the 
fuel  cell  does  not  justify  increasing  the  flow  rate.  The  power  output 
of  the  compressor  is  given  by  [7] 


Current  =  50  mA  cm  2 


5  10  IS  20  25  30x10* 

Oxygen  Flow  Rate  /  mol  s’ 


where  y  is  the  heat  capacity  at  constant  pressure  divided  by  that  at 
constant  volume  (y  =  1.4),  is  the  inlet  oxygen  pressure  in  the 
cathode  flow  channel,  ECOmp  is  the  air  compressor  efficiency 
(assumed  to  be  0.8)  and  the  coefficient  of  371  has  units  of  J  K-1  m~3. 

The  blower  power  output  can  be  found  using  the  specific  fan 
power  (SFP)  as  follows: 

Powerblower  =  SFP*JV®put.  (36) 

The  lowest  SFP  values  for  small  commercial  fans  operating 
under  normal  conditions  is  approximately  500  W  s  m  3  [29],  so  we 
choose  to  use  that  value  through  this  analysis.  Note  the  need  to 
convert  the  oxygen  flow  rate  from  mol  s-1  to  m3  s_1  through 
multiplication  by  the  molar  volume.  Although  the  compressor  uses 
much  more  power  than  the  blower,  the  performance  gain  achieved 
by  raising  the  cathode  input  pressure  is  substantial  and  the  need 
to  maintain  a  high  oxygen  concentration  in  the  presence  of  meth¬ 
anol  crossover  is  imperative.  Therefore,  we  chose  to  include 
a  compressor  throughout  this  analysis.  Subtracting  these  parasitic 
losses  from  the  power  output  of  the  fuel  cells  provides  us  with  an 
expression  for  the  net  power  output  of  the  system 


Powersystem  =  Vcell/Ceii(#cells)  -  Powercomp  -  Powerblower 


-  Powerpump. 


(37) 


Because  the  analysis  below  is  only  concerned  with  the  oxygen 
flow  rate,  the  parasitic  losses  from  the  methanol  pump  Powerpump 
are  ignored  in  the  calculation  of  the  system  power  output. 

Fig.  8  above  shows  the  power  output  of  a  DMFC  stack  containing 
100  5  cm2  cells  as  a  function  of  oxygen  flow  rate.  As  expected,  the 


Fig.  8.  Net  power  output  of  a  100  cell  stack  as  a  function  of  oxygen  flow  rate  with  0.05 
and  0.3  A  cnr2  currents  using  Nafion  115  and  C-10  blend  membranes.  All  other 
parameters  and  operating  conditions  are  identical  to  what  is  used  in  Fig.  3. 


80  |im  thick  blend  membranes  outperform  the  Nafion  115 
membranes,  and  the  1  M  MeOH  feed  stream  is  always  better  for 
performance  than  2  M  MeOH  due  to  zero-order  methanol  kinetics. 
Similar  to  what  was  done  in  Fig.  7,  the  desired  current  is  lowered  at 
low  flow  rates  in  order  to  keep  the  oxygen  pressure  above  zero 
across  the  entire  cathode.  However,  in  some  cases  at  50  mA  cm-2 
the  current  cannot  be  lowered  any  further,  and  so  the  fuel  cell 
cannot  run  at  all  below  those  cutoff  flow  rates.  The  2  M  Nafion  115 
case  in  particular  shows  a  much  higher  cutoff  flow  rate  than  the 
blend  membranes.  This  low  current  scenario  is  where  the  blend 
membranes  have  the  greatest  advantage  over  the  Nafion 
membranes. 

In  all  of  these  cases,  the  power  demand  of  the  compressor  and 
the  blower  is  a  substantial  fraction  of  the  power  output  of  the  fuel 
cell  stack.  At  the  optimum  flow  rate  in  the  300  mA  cm-2  scenario, 
the  power  output  of  the  100  cell  stack  is  about  40  W,  and  the 
power  consumption  of  the  compressor  and  the  blower  is  about 
4  W,  or  10%  of  the  total  power  output.  That  fraction  only  gets 
larger  as  the  flow  rate  increases  because  the  performance  gains 
become  smaller  and  smaller  while  the  power  consumption  of  the 
system  components  continues  to  increase  linearly.  The  oxygen 
reduction  reaction  is  first-order,  so  there  is  a  noticeable  perfor¬ 
mance  penalty  to  be  paid  for  operating  close  to  the  cutoff  flow 
rate.  However,  until  the  flow  rate  is  within  about  2  x  10-6  mol  s-1 
of  the  cutoff  flow  rate,  that  penalty  is  more  than  compensated  for 
by  the  lower  power  consumption  of  the  system  components. 
Therefore,  as  with  the  methanol  flow  rate,  it  is  important  to  keep 
the  oxygen  flow  rate  as  low  as  possible  while  not  allowing  the 
oxygen  pressure  to  fall  close  to  zero  such  that  substantial  perfor¬ 
mance  losses  are  incurred. 
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5.  Conclusion 

In  this  work,  a  simple  engineering  model  is  developed  that 
directly  relates  the  materials  properties  of  the  direct  methanol  fuel 
cell  —  including  diffusion  coefficients,  membrane  thickness,  and 
proton  conductivity  —  and  operating  parameters  —  most  notably 
the  methanol  and  oxygen  flow  rates  —  to  the  stack  performance  and 
the  size  of  system  components  such  as  the  methanol  storage  tank 
and  the  air  blower.  The  model  is  created  by  deriving  analytic 
solutions  for  the  methanol  and  oxygen  concentrations  across  the 
thickness  of  the  fuel  cell  and  by  applying  Tafel  kinetics  to  determine 
the  average  cell  current  density  for  a  fixed  cell  voltage.  The  benefits 
of  the  model  include  a  run  time  of  less  than  a  minute  for  a  single 
current-voltage  curve  in  a  MATLAB  simulation  environment  and 
the  ability  to  handle  a  wide  range  of  parameter  spaces. 

The  model  was  validated  using  experimental  current— voltage 
data  taken  at  65  °C  and  90  °C  with  air  oxygen  cathode  feeds.  It 
was  then  applied  to  forecast  the  benefits  of  using  new  low- 
crossover  membranes  and  was  compared  to  published  data  from 
Zhu  et  al.  Finally,  the  model  was  used  in  a  case  study  to  illustrate  the 
effect  of  changing  the  methanol  and  oxygen  flow  rates  on  fuel  cell 
performance  and  the  net  power  output  of  the  fuel  cell  stack.  The 
tradeoff  that  exists  between  stack  performance,  the  methanol 
storage  tank,  and  the  power  consumption  of  the  air  blower,  air 
compressor,  and  other  components  dictates  that  there  exists  an 
optimum  set  of  operating  conditions  for  any  portable  DMFC  with 
size  constraints. 
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Glossary 


Roman  Symbols 

A  superficial  area  of  the  anode,  cm2 

Achan  cross-sectional  area  of  the  flow  channel,  cm2 
Canode  methanol  concentration  at  the  anode  catalyst  layer, 
mol  cm-3 

Cfeed  methanol  concentration  in  the  anode  feed  stream,  mol  cm  3 

cfeed  methanol  feed  concentration  at  channel  inlet  (z  =  0), 
mol  cm-3 

Ch2o  concentration  of  liquid  water  in  the  membrane,  mol  cm”3 

Dj  _  j  diffusion  coefficient  of  component  i  in  component,  cm2  s  1 
Dflfj  effective  diffusion  coefficient  of  component  i  in 
component  j,  cm2  s”1 

Dmem  diffusion  coefficient  of  methanol  in  the  membrane,  cm2  s”1 

Dq2  ci  Deff  of  oxygen  in  the  cathode  catalyst  layer,  cm2  s”1 
df,  hydraulic  diameter  of  the  flow  channel,  cm 

Ecomp  air  compressor  efficiency,  dimensionless 

/  friction  coefficient  of  air  in  the  flow  channel, 

dimensionless 

h  size  of  each  iteration  step  in  the  flow  channels,  cm 

i  current  through  the  external  circuit,  A  cm-2 

ix  crossover  current,  A  cm-2 

ig  anodic  reference  exchange  current  density,  A  cm-2 

cathodic  reference  exchange  current  density,  A  cm  2 
iiim,a  limiting  current  in  the  anode  backing  layer,  A  cm-2 
'Em  a  backing  layer  limiting  current  at  anode  channel  inlet, 

C  cm-2  s”1 

k  cross-sectional  area  to  flow  in  cathode  flow  channel,  cm2 

f-back,a  thickness  of  the  anode  backing  layer,  cm 
f-back,c  thickness  of  the  cathode  backing  layer,  cm 


F-mem  thickness  of  the  membrane,  cm 
Id  thickness  of  the  catalyst  layers,  cm 
Ichan  length  of  the  flow  channels,  cm 
n  total  number  of  iterations  down  the  flow  channel, 

dimensionless 

Nj  flow  rate  of  component  i  through  region  j,  mol  s”1 

Nif  flow  rate  of  component  i  in  the  cathode  flow  channel, 

mol  s”1 

NMeOH  methanol  flow  rate  into  the  anode  flow  channel  inlet, 
mol  s”1 

Ninput 

oxygen  flow  rate  into  the  cathode  flow  channel  inlet, 
mol  s”1 

P,  absolute  pressure  of  component  i,  N  m  2 

Pi/  pressure  of  component  i  in  the  bulk  flow  channel,  N  m”2 

Po2  s  oxygen  pressure  at  the  cathode  backing  layer  surface, 

N  nr2 

Po2,c  oxygen  pressure  at  the  cathode  catalyst  layer  surface, 

N  m”2 

Po2.ref  reference  oxygen  pressure,  N  m  2 

Pq2  j  input  oxygen  pressure  in  the  cathode  flow  channel,  N  m”2 

Ph2o  sat  saturated  vapor  pressure  of  water,  N  m”2 

pbatk 

oxygen  pressure  in  the  cathode  backing  layer,  N  m  2 
Pq'11v  oxygen  pressure  in  convective  flow  channel  regions, 

N  m”2 

Ptotai  total  pressure  in  the  cathode  backing  layer  N  m”2 
PowerComp  power  consumed  by  the  air  compressor,  W 
Powerbiower  power  consumed  by  the  air  blower,  W 
Re  Reynolds  number  of  cathode  air  flow,  dimensionless 
Si  stoichiometric  ratio  of  methanol  or  oxygen, 

dimensionless 

Shp  Sherwood  number  for  a  square  channel,  dimensionless 

SFP  specific  fan  power,  W  s  m  3 

t  time  of  fuel  cell  operation,  s 

if  standard  cell  potential,  V 

kanode  anodic  activation  overpotential,  V 

^cathode  cathodic  activation  overpotential,  V 

Vceii  total  cell  potential,  V 

Vfeed  molar  volume  of  the  methanol  and  water  feed  stream, 

cm3  mol-1 

VMeOH  velocity  of  methanol  in  the  anode  flow  channel,  cm  s”1 

Vair /  velocity  of  air  in  the  cathode  flow  channel,  cm  s”1 

vo2/  velocity  of  oxygen  in  the  cathode  flow  channel,  cm  s”1 

wsep  flow  channel  separation  width,  cm 

yi  mole  fraction  of  component  i,  dimensionless 

z  distance  down  the  flow  channel,  cm 

z'  normalized  distance  down  the  flow  channel, 

dimensionless 


Creek  Symbols 

aa  anodic  transfer  coefficient,  dimensionless 

ac  cathodic  transfer  coefficient,  dimensionless 

ea  porosity  of  the  anode  backing  layer,  dimensionless 
ec  porosity  of  the  cathode  backing  layer,  dimensionless 

£q\  porosity  of  the  catalyst  layers,  dimensionless 
kmem  proton  conductivity  of  the  membrane,  S  cm-1 
Kcia  proton  conductivity  of  the  anode  catalyst  layer,  S  cm  1 

k dc  proton  conductivity  of  the  cathode  catalyst  layer,  S  cm”1 

Hi  dynamic  viscosity  of  component  i,  Pa  s”1 
t]  effectiveness  factor  of  cathode  catalyst  layer, 

dimensionless 

^Nafion  Nafion  volume  fraction  in  the  catalyst  layer, 
dimensionless 

tj>i  sum  of  dimensionless  ratios  (see  pp.  153—155  in  Ref.  [23]), 

dimensionless 
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\p  ratio  of  ifjm  to  the  current  through  the  external  circuit, 

dimensionless 

Ph2o  density  of  liquid  water,  g  cm  3 
8  Thiele  modulus,  dimensionless 

v  ratio  of  limiting  current  in  backing  layer  and  membrane, 

dimensionless 

?h2o  electroosmotic  drag  coefficient  of  water,  cm3  mol-1 
?MeOH  electroosmotic  drag  coefficient  of  methanol,  cm3  mor1 
?MeOH  electroosmotic  drag  factor  of  methanol,  dimensionless 
7  heat  capacity  ratio  (7p/yv)  dimensionless 
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